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Abstract —Inverse determinant sums appear naturally as a 
tool for analyzing performance of space-time codes in Rayleigh 
fading channels. This work will analyze the growth of inverse 
determinant sums of a family of quasi-orthogonal codes and 
will show that the growths are in logarithmic class. This is 
considerably lower than that of comparable number field codes. 

I. Introduction 

In Q inverse determinant sums were proposed as a tool 
to analyze the performance of algebraic space-time codes 
for MIMO fading channels. These sums can be seen as a 
generalization of the theta series for the Gaussian channel. 
They arise naturally from the union bound on the pairwise 
error probability for spherical constellations, but also in the 
analysis of fading wiretap channels f6j. 

In 0 the authors analyzed the growth of the inverse deter¬ 
minant sums of diagonal number field codes and of most well 
known division algebra codes. In this work we are going to 
extend the analysis to a large class of quasi-orthogonal codes. 
Our work will reveal that the growth of inverse determinant 
sums of the analyzed codes is considerably smaller than that of 
the corresponding diagonal number field codes. This difference 
suggest that asymptotically, with growing constellation, quasi- 
orthogonal codes are considerably better than number field 
codes. This difference can not be captured in the framework 
of diversity-multiplexing gain tradeoff. 

For related work, we refer the reader to (T) and fD- 

Remark 1.1: This is a corrected and extended version of 
0. The main correction is to Theorem 13.41 which we quoted 
carelessly in the previous version. This correction does affect 
all the formulas that are derived from it. However, it does not 
change any formulas related to complex constellations. In the 
case of totally real fields the change sometimes doubles the 
number of units. The other correction is to the definition of 
the normalized inverse determinant sum Q. This correction 
does not change any results, but is corrected for consistency. 

The extensions in this paper are Propositions 13.81 and 13.61 
These results follow from improving on the estimate for 
the truncated Dedekind zeta function done in 0- The new 
estimate is Lemma [3731 The proof of this result is due to Tom 
Meurman and the authors are grateful to him. 

Remark 1.2: We note that after the publication of [[2] the 
authors in ITOl extended Proposition 13.51 to ideals and proved 


that our upper bound is tight in the case of principal ideal 
domains. 


II. Inverse determinant sum 

We begin by providing basic definitions concerning matrix 
lattices and spherical constellations, that are needed in the 
sequel. 

A. Matrix lattices and spherically shaped coding schemes 

Definition 2.1: A space-time lattice code C C M„(C) has 
the form 

ZHi © © • • • © 

where the matrices B \,..., B k are linearly independent over 
M, i. e. , form a lattice basis, and k is called the rank or the 
dimension of the lattice. 

Definition 2.2: If the minimum determinant of the lattice 
L C M„(C) is non-zero, i.e. it satisfies 

inf |det(X)| > 0, 

o^A'ei 

we say that the code has a non-vanishing determinant (NVD). 

We now consider a coding scheme based on a fc-dimensional 
lattice L inside M r: (C). For a given positive real number M 
we define the finite code 

L(M) = {a | a G L, ||a|| F < M}, 

where ||a|| F refers to the Frobenius norm. In the following we 
will also use the notation 


B{M) = {a | a £ M„(C), ||a|| F < M}, 


for the sphere with radius M. 

Let L C M n (C) be a fc-dimensional lattice. For any fixed 
to G Z+ we define 


S?(M) := 


E 

X£L(M)\{ 0 } 


1 

| det(X)| m ’ 


Our main goal is to study the growth of this sum as M 
increases. Note, however, that in order to have a fair com¬ 
parison between two different space-time codes, these should 
be normalized to have the same average energy. Namely, the 
volume Vol(fc) of the fundamental parallelotope 


V(L) = {aiBi + a 2 B 2 + .. ■ + a k B k | cu G [0,1) Vi} 



should be normalized to 1. The normalized version of the 
inverse determinant sums problem is then to consider the 
growth of the sum S™(M) = S™(M) over the lattice L = 
Vol {L)~ 1 / k L. Since L(M) = Vo\{L)- 1 / k L(M Vol(X) 1 / fe ), 
we have 

S™(M) = Vo\{L) mn/k S?(MVol{L) 1/k ). (1) 
Remark 2.1: In El we erroneously stated that 

S™(M) = Vol (L) mn ' k S^{M). 

The correct form is ©, which was originally stated in The 
difference is that in the correct expression the radius of the 
ball also increases. However, this difference does not affect 
the results in this paper. In every case in this paper the growth 
of the inverse determinant is logarithmic and therefore the 
constant term Vol (L) l / k will have no effect on the dominating 
term of the sum. 


B. Inverse determinant sums and error performance of space- 
time lattice codes 

Let us now consider the slow Rayleigh fading MIMO 
channel with n transmit and n r receive antennas. The channel 
equation can then be written as 

Y = HX + N , 


where H and N are respectively the channel and noise 
matrices. We suppose that the transmitted codeword X belongs 
to a finite code L(M) C M n ( C) carved from a fc-dimensional 
NVD lattice L as defined previously. In terms of pairwise error 
probability, we have for I/V 


P(X -s- X') < 


1 

| det(V — X')\ 2nr ’ 


and the corresponding upper bound on overall error probability 


Pe < 


E 

X£L,0<\\X\\ f <2M 


1 

|det(V)| 2Ti V 


Our main goal is now to study the growth of the sum 
S™{M) as M increases. In particular, we want to find, if 
possible, a function f(M) such that 


S?(M) ~ f(M). 


III. Inverse determinant sums of algebraic 

NUMBER FIELD CODES 

In this section we will give and review some results concern¬ 
ing inverse determinant sums of diagonal number field codes. 
These results will play an important role in our analysis of 
quasi-orthogonal codes. The proofs are mostly analogous to 
those given in Q and we will skip them. However, the case 
where the receiver has a single antenna will improve on the 
original one given in (7j. Unlike in the rest of the paper we 
will state the results in the normalized from S™(M) following 
the general normalization given in 0. 

We need first some preliminary results. Let us suppose we 
have an algebraic number field K. We will denote with Ok 


the ring of algebraic integers in K. If B is an integral ideal 
of K we will use the notation N(B) = [Ok ■ B]. 

Let Iq k be the set of nonzero ideals of the ring Ok- The 
Dedekind zeta function of the number field I\ is 

C * (S)= E 7 ’ (2) 

Aeio K y ’ 

where s is a complex number with Vt(s) > 1. For a truncated 
version of this sum we will use notation 

C k{s,M):= Y, 

A£I 0 k ,N(A)<M y ’ 


In the following we will need an estimate for (k(s,M). We 
will begin with a classical result. Let us suppose we have an 
algebraic number field K with signature 

Theorem 3.1: The function £ k{s) converges for all s, 
5R(s) > 1 and has a simple pole at 1 with residue 


lim(s — l)Of(s) 

S-5-1 


2^(2 x) r *h K R K 

uVWk)\ 


where R k is the regulator of the number field K , d(K) the 
discriminant and ui the number of roots of unity in K. 

In the following we will use the shorthand 

2 ri (27T Y 2 Rk 
OIK = - , 

uVmT\ 

In order to estimate the truncated Dedekind zeta function at 
the point 1 we need an estimate for the number of ideals with 
bounded norm. Let us denote with N(K. M) the number of 
integral ideals of K with norm less than or equal to M. We 
then have the following. 

Proposition 3.2 : 0 Theorem 5] Given a degree n number 
field K there exist positive constants c and a, independent of 
M, such that 


\N(K, M) - a K h K M\ < cM 1 ^. 


We are now ready to estimate the truncated Dedekind zeta 
function. 

Lemma 3.3: We have the following 

E nTaS = h K a K logM+ 0(1). 

AgIo k ,N(A)<M y ' 


Proof: Let us first write 


E 

N(A)<M 


1 

N(A) 


E l E i 

n<M N(A)=n 


E 

n<M 


z{n) 

n 


Manipulating the previous expression, we can write it in the 
form 


S + T= Y 

n<M 


z(n) - h K oiK 
n 


+ E 

n<M 


h-K^X 

n 


(3) 


Let us now analyze the first term S by using partial summation 
U2 A. 21] 

1 r M i 

S' = — E ( Z ( n ) - h KOtx) + 72 E(~( n ) “ h KOi K )dt. 

n<M n<t 















According to Proposition 13.21 we have that 

^2(z{n) - h K a K ) = ^2 ~( n ) “ h KUKt + 0(1) = 0(f 1_a ). 


n<£ 

It follows that 

r M 


n<t 


/ M i pM 

— '22(z(n)-h K a K )dt= 0(f _1_a )df. 

n<t •> 1 

Putting now all together we have that 
i r M 

S = —0(M 1 ~ a ) + / 0(t~ l ~ a )dt = 0(M~ a ) + 0(1). 
M J i 

Estimating J2 n <M hj *7T L = h K a K log M + 0(1), gives 


us 


the final result. □ 

Remark 3.1: This results improves on the very crude esti¬ 
mate 0 

1 


E 


A£Io k ,N(A)<M 

where d is some constant. 


N(A) 


< d(logM)", 


A. Inverse determinant sums of real diagonal number field 
codes 

Let A' be a totally real number field of degree n and let 
{cti, - - - ,cr„} be the Q-embeddings from K to R. We then 
have the canonical embedding ip : K n- M„(R) defined by 

f(x) = diag(ui(a;),..., a n (x)). 

It is a well known result that ^{Ok) is an /(-dimensional 
NVD lattice in Let us now consider the corresponding 

inverse determinant sum. The main role in the analysis is 
played by the following unit group density result. 

Theorem 3.4 (K8f): Let us suppose that [A' : Q] = n, we 
then have that 


where c K = ^-lf ■ 

Proof: Lollowing the proof of Proposition 4.4 and Section 
4.C in 0 we have an estimate 

M n 

S™ ( g k) {M) < Gf(l, —^){N K (logM)™- 1 +0((\ogM) n ~ 2 )), 

and the end result now follows from plugging the result of 
Lemma [3731 and Theorem [+4] to this estimate. □ 

Lor principal ideal domains the upper bound is likely tight. 

B. Inverse determinant sums of complex diagonal number field 
codes 

Let AT/Q be a totally complex extension of degree 2 n and 
{<Ti,..., a n } be a set of Q-embeddings, such that we have 
chosen one from each complex conjugate pair. Then we can 
define a relative canonical embedding of AT into M n (C) by 

ip(x) = diag(ai(x),... ,a n (x)). 

The ring of algebraic integers Ok has a Z-basis W = 
{wi,... ,W 2 n} an d f’iW) is a Z-basis for the full lattice 
iP{O k ) in M„(C). 

Proposition 3. 7: 0 Section 4.C] Let K be a totally com¬ 
plex algebraic number field of degree 2 n. If n r > 1, we have 
that 

^ k) (M) < NkCkM (log M) n ~ l + 0((log M) n ~ 2 ) 

and 

N k (logM)"" 1 + 0((logM)"- 2 ) < S^ k) (M), 
where N K = §^(2 

Proposition 3.8: Let AT be a totally complex algebraic 
number field of degree 2 n. We then have 

S* ( o*)(M) < ck (log M) n + O(aogM)"- 1 ) (4) 

and 


W(0* K )nB(M)\ = N K (log My - 1 +0((logM) n ~ 2 ), 
where N K = 

Here R is the regulator of the number field AT, ui the number 
of roots of unity in K. 

We first have a real analogue to the results in 0. 
Proposition 3.5: Let us suppose that K is a totally real 
number field with [K : Q] = n and that m > 1. Then 


S$ ( o K )(M) < N K ( K (m) (log M) r 

and 


0((logM) r 


N k (logM)"" 1 + 0((log M) n ~ 2 ) < S^ 0 k) (M), 

where N K = 0Sf(VWKT\) m - 

In the case of single receive antenna we now have the 
following. 

Proposition 3.6: Let us suppose that K is a totally real 
number field with [AT : Q] = n and that m = 1. Then 

S^ { o k) (M) < CK (log Mf + 0((logM)™ -1 ), 


N k (log M) n ~ x + 0{(logM) n ~ 2 ) < S 2 {0 k) (M), (5) 

where ^f and N K = gg^J(2 

Proof: Lollowing the proof of Proposition 4.4 in 0 we 
have an estimate 

M 2n 

S^ 0k) (M) < C x(l,^ r )(iVi f (logM)"- 1 +0((logM)"- 2 ), 

and the end result now follows from plugging the result of 
Lemma [3+1 and Theorem 13.41 to this estimate □ 

Again for principal ideal domains the upper bound is likely 
tight. 

IV. Quasi-orthogonal codes from division 

ALGEBRAS 

In the following we are considering the Alamouti-like 
multiblock codes from 0. With respect to their complexity 
and other properties, all of the codes of this type are quasi- 
orthogonal. It is even possible to prove that many of the fully 
diverse quasi-orthogonal codes in the literature are unitarily 
equivalent to these multi-block codes. In the following we 















will use several results and concepts from the theory of central 
simple algebras. We refer the reader to fill for an introduction 
to this theory. 

Let us consider the held E = KF that is a compositum 
of a complex quadratic held F and a totally real Galois 
extension K /Q of degree k. We suppose that K fl F = Q, 
Gal(F/Q) =< a > and Gal{K/ Q) = {n, r 2 ,..., 77 }. Here 
a is simply the complex conjugation. We can then write that 
Gal(FK/Q) = Gal(K/ Q)<g> < cr >. 

Let us now consider a cyclic division algebra 

V = ( E/K , cr, 7 ) = i? ® uE , 


where u £ V is an auxiliary generating element subject to the 
relations xu = ux* for all x £ E and u 2 = 7 £ Ok* where 
()* is the complex conjugation. We can consider T> as a right 
vector space over E and every element a = x\ + 11 x 2 £ V 
maps to 


</>(a) 


Xl 

1X* 2 



This mapping can then be extended into a multi-block 
representation <j>: V i-> M 2 fe(C). 


= diag(ri(</>(a)),r 2 (0(a))..., r fe (</>(a))). (6) 

nple 4.1: In the case where k = 2 each element a £ V 

x 2 0 0 \ 

xl 0 0 

0 t(x 1 ) t(x 2 ) 

0 t( r yx 2 )* t(xi)*J 

In order to build a space-time lattice code from the division 
algebra D we will need the following definition. 

Definition 4.1: Let Ok be the ring of integers of I\. An 
OK-order A in V is a subring of V, having the same identity 
element as V, and such that A is a hnitely generated module 
over Ok and generates V as a linear space over I\. 

Let us suppose that A is an Off-order in V. We call <p(A) 
an order code. In the rest of this paper, we suppose that the 
division algebras under consideration are of the previous type. 

Lemma 4.1: If A is an Off-order in V 

| det((j)(x))\ = a/[A : xA], (7) 




gets mapped as 




( Xi 

0 

0 


where a: is a non-zero element of A. 

Lemma 4.2: Let us suppose that A is a Ok- order of a 
division algebra V with center K of degree k and that <p is 
a multi-block representation. Then the order code (f>( A) is a 
4fc-dimensional lattice in the space M 2 fc(C) and 

dot rnin (t/ ? (A)) = 1. 


Let V be an index-n A'-central division algebra and A a 
Ok- order in V. The (right) Hey zeta function of the order A 
is 


Ca(s) = E 

teix 


1 

[A : J] s ’ 


where Ji(s) > 1 and Ia is the set of right ideals of A. When 
5ft(s) > 1 , this series is converging. 

The unit group A* of an order A consists of elements x £ A 
such that there exists a y £ A with xy = 1^. Another way to 
define this set is A* = {x £ A | | defeat) | = 1}. 


A. Inverse determinant sums of quasi-orthogonal codes 
Let us suppose that K, V and A are as in the previous 
section and that [K : Q] = k. We then have that <j>(A) is a 
4fc-dimensional NVD lattice in M 2 k (C) and we can consider 
the growth of the sum 


E 


1 

| det-!/>(x)| 2rir 


_ q2n r 

- ^(A) 


(M). 


Just as in (71 the previous sum can be analyzed further into 


q2n r 

*v(A) 


(M) 


x ^ \il>(xA*)r\B{M)\ 

I det(^>(a;))| 2nr 


( 8 ) 


where X(M) is some collection of elements x £ A such that 
||V>(ar)||p, < M, each generating a different right ideal. 


B. Uniform upper and lower bounds for ^(xA*) H B(M)\ 
The key element in the analysis of \i/j(xA*) C\B(M) \ is the 
following. 

Lemma 4.3 (Eichler): The unit group A* has a subgroup 


0* K = {x | x £ A*, x £ O k }, 

and we have [A* : O < oo. 

Let j = [A* : Ok]* By choosing a set {ai,..., a 7 } of coset 
leaders of 0* K in A*, we have that 

3 

\ip{xA*) n B(M) | < ^ | fi(xa t O* K ) n B(M)\. ( 9 ) 

i=l 

In order to give an uniform upper bound for \i/j(xA*) n 
B(M)\, it is now enough to give a uniform upper bound for 
\fi>{xaiO* K ) C\B(M)\. Before stating our main results we need 
few lemmas. We will skip the proofs of some of them. 

Lemma 4.4: Let us suppose that A is a diagonal matrix in 
M n (C) with | det A\ > 1. We then have that 

\Af{0* K )GB{M)\ < \fi{0* K ) (T B(cM)\, 

where c is a fixed constant, independent of A and M. 

Lemma 4.5: Let us suppose that x and y are elements in 
Okf* we then have that 


\fj(xMO* K )nB(M)\ < \^(p* K ) n b{ c m)\, 


and 


Muy)i>(0* K )nB(M)\ < \i/j( 0* k ) n B(cM)\, 

where c is a real constant independent of x, y and M. 

Proof: The first result is simply Lemma 14.41 and the 
second follows as ip(u) is a fixed matrix. □ 

Lemma 4.6: Let us suppose that x and y are elements in 
E. We then have that 


\\ft{x) +tp(uy)\\ 2 F = \\ip(x)\\ 2 F + \\ip{uy)\\ 2 F . 








Proof: By an elementary calculation we see that < 
ip(x),ip(uy) >= 0 and the claim follows. □ 

Proposition 4.7: Let us suppose that x £ A, we then have 
that 


n B(M )| < mo* K ) n B(cM) I, 

where c is a constant independent of M and x. 

Proof: Let us suppose first that x = x\ + ux 2 , where 
Xi £ O e and where u 2 £ Ok- According to Lemma [4.61 we 
have that 

\\xl){x)f){y)\\ 2 = ||0(a;i)^(t/)|| 2 + \\4’(ux 2 )ip(y)\\ 2 , 
for any y £ Oe- 

Therefore if ip(x)ip(y) £ B(M), then also 

t/)(xi )ip(y) £ B(M) and f(uX 2 )ip{y) £ B(M). 

It follows that we can upper bound | ip(x)ip(0* K ) H B(M)\ 
with 

max{| fj( Xl )fj(0* K ) n B(M )|, | i>(ux 2 )f>(0* K ) n B(M )|}. 

According to lemma 1431 we then have that 

MxMO* k ) n B{M )I < I 4>(o* k ) n B(cM)\. 

Let us now suppose that A is a general order in T>. As ip(A) 
is finitely generated as an additive group in M n (E), we can 
choose an integer d such that dip (A) C iP(Oe) + iP{uOe )• 
The result now follows from the previous consideration. □ 
Proposition 4.8: Using the previous notation we have 

\ip(xA*) n B(M)\ < [A* : O k ] ■ 

log M k ~ x k _ ^ + 0(log M k ~ 2 ). 

Proof: Let j = [A* : Ok] - By choosing a set {ai,..., a,j} 
of coset leaders of 0* K in A*, we then have that 

j 

| ip{xA*) n B(M )| < \^{xaiO* K ) n B(M)\. 

i=l 

According to Proposition 14.71 we then have that 

\iP(xA*)nB(M)\ < [A* :0 K ]MO* K )nB(cM)\. 

Applying Theorem l3.4l to this equation, we get the final result. 

□ 


C. Upper and lower bounds for inverse determinant sums of 
quasi-orthogonal codes 

Proposition 4.9: Let us suppose that [K : Q] = k and set 
n = 2k. We then have that ip( A) is a 2n-dimensional lattice 
in M n (C) and 


. ,( n \ n /2 —1 
t /2 —1 U \ 2 > _ 




< CaK)[A* : Ok] logM^ 


R{^)\ 


+ 0(\ogM n/2 ~ 2 ), 


where n r > 1 and R and uj are the regulator and the number 
of roots of unity in the center K. 

Proof: As previously mentioned, we can imitate Q to 


get 



E 

xex(M) 


\ip{xA*)nB{M)\ 
| det(^(a;))| 2nr 


( 10 ) 


According to Lemma |4~TI we have that |det(t)’(a;))| 2n ' r = [A : 
xA] nr . Now 

,Jiv, iSwwir £ ^ £ CaK) ' 

Applying this inequality with Proposition l4.8l to ( fTTTb now gives 
us the final result. □ 


V. Quasi-orthogonal codes are better than 

DIAGONAL NUMBER FIELD CODES 

Let us now suppose we have an n x n r -MIMO channel, 
(for simplicity we assume n r > 1). For the existence of quasi- 
orthogonal code we also have to assume that 2 | n. Let us now 
compare the growth of determinant sums of quasi-orthogonal 
and comparable diagonal number field codes in this n x n r - 
MIMO channel. 

In order to build a quasi-orthogonal code ip (A) in Af„(C) 
the center K of the algebra T> must be an n/2-dimensional 
totally real number field. For a number field code ip{0]f) C 
M„(C), the field L must be an n-dimensional extension of 
some complex quadratic field F. 

As we earlier saw, we have that 

L ^TF = w<A ‘ )nB(M)l) 

Xeif>(A)(M) 1 v n 

and 

]ip{A*) n B(M )I = 6 »(| ip{0* K ) n B(M) I) = 0(log AT/ 2 " 1 ). 


Therefore 


E 


1 


|det(AT)| 2 ^ 


= 0(logAf n/2 " 1 ). 


Xeip(A){M) 

On the other hand for the number field code we have that 
1 


E 


xeip(o L )(M) 


|det(X)| 2 "r 


= e(\ip(o* L )nB(M)\) 


= 0(log AI 71-1 ). 


Here the last result follows from 0 Theorem 2], 

We can now see that the growth of the inverse determinant 
sum for the quasi-orthogonal code is considerably lower than 
that of the number field code. This is due to the fact that the 
unit group of the order A is essentially that of a low degree real 
number field. We note that this difference can not be captured 
in the context of DMT as both of these codes have the same 
DMT curve. 
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